The Anatomy of Star Formation in NGC 300 
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ABSTRACT 

The Spitzer Space Telescope was used to study the mid-infrared to far-infrared 
properties of NGC 300, and to compare dust emission to Ha to elucidate the 
heating of the ISM and the star formation cycle at scales < 100 pc. The new 
data allow us to discern clear differences in the spatial distribution of 8 //m dust 
emission with respect to 24 /im dust and to H II regions traced by the Ha light. 
The 8 /im emission highlights the rims of H II regions, and the 24 /im emission 
is more strongly peaked in star forming regions than 8 /xm. We confirm the 
existence and approximate amplitude of interstellar dust emission at 4.5 fim, 
detected statistically in Infrared Space Observatory (ISO) data, and conclude it 
arises in star forming regions. When averaging over regions larger than ~ 1 kpc, 
the ratio of Ha to Aromatic Peature emission in NGC 300 is consistent with the 
values observed in disks of spiral galaxies. The mid-to-far-infrared spectral energy 
distribution of dust emission is generally consistent with pre-Spitzer models. 

Subject headings: galaxies: ISM — infrared: galaxies — individual: galaxy 
(NGC 300) 

1. Introduction 

NGC 300 is a SA(s)d galaxy in the Sculptor Group of galaxies, at a distance of about 
2.1 Mpc (Preedman et al. 1992), viewed at an inclination of about 50° (Puche et al. 1990). 
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Its total luminosities in the blue band (3855-4985A) and in the far-infrared (42.5-122.5 //m, 
estimated from Spitzer data) are ~ 3.3 x 10^ Lq i,oi and 2.2 x 10^ i,oi respectively. Its 
Lfir/Lb ratio is thus very close to the average ratio for the blue magnitude-limited sample 
studied by Thuan & Sauvage (1992). Due to its large angular extent and low surface bright- 
ness, no reliable total radio continuum flux measurement exists (only discrete sources are 
detected). NGC 300 has a striking appearance in the visible and in Ha due to several H II 
regions with nearly circular shapes and various degrees of fiUing-in (Deharveng et al. 1988; 
Hoopes et al. 1996; see Figure 1). Its large HI envelope (Puche et al. 1990) extends well 
beyond the visible image. Pannuti et al. (2000) reported a total of 44 SNR candidates, evi- 
dence that the current star formation activity manifested as H II regions has been on-going 
for tens of millions of years (see also Butler et al. 2004). Because of its proximity, NGC 300 
allows the Spitzer Space Telescope (Werner et al. 2004) to compare dust emission to other 
ISM components and discern the interplay between the ISM and the star formation cycle at 
scales < 100 pc. 



2. Observations and Data Reduction 

The IRAC (Fazio et al. 2004) observations of NGC 300 are 12-second frames and map 
the galaxy in approximately half array spacings, yielding a total time per sky position of 48 s. 
They were reduced with the standard Spitzer Science Center data reduction pipeline (version 
9.5). Due to the readout of bright point sources, the signal was reduced (pulled-down) in some 
columns of the array^; a correction measured from the map was applied. Persistent images 
left by bright sources were found by median-combining all the dithered positions together 
and identifying any remaining sources. These objects were removed from individual frames. 
The data were then combined into a mosaic using a cosmic ray rejection and a background 
matching applied between overlapping fields of view. The relative photometric uncertainty 
is of the order of 5%, and the absolute uncertainty is 10%; uncertainty tied to the angular 
sizes of source and measuring aperture contributes up to 15%. 

Images of NGC 300 at 24, 70, and 160 //m were obtained with the MIPS Instrument 
(Rieke et al. 2004) in the scan-map mode. The final mosaics have a total exposure time 
of approximately 160, 80, and 16 s per point at 24, 70, and 160 /im. The MIPS images 
were reduced using the MIPS Instrument Team Data Analysis Tool (Gordon et al. 2004) as 
described by Engclbracht et al. (2004). The uncertainties on the final absolute calibrations 
arc estimated at 10%, 20%, and 20% for the 24, 70, and 160 /xm data, respectively. The 
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70 and 160 //m images exhibit linear streaks along the scan direction which are residual 
instrumental artifacts due to the time dependent responsivity of the Ge detectors, and affect 
the photometry in large apertures. 

The Ha (6563A) + [N II] (6583 and 6548A) map was derived from images posted in 
NED by Larsen & Richtler (1999). The narrowband image containing Ha emission and the 
R-band image were aligned, and rescalcd to subtract the stellar continuum, deriving the 
scaling factor from aperture photometry on 32 bright stars. Residuals from saturated or 
improperly subtracted stars were masked out, and the final image was regridded and rotated 
to match the IRAC reference frame. 

The Spitzer and Ha+[N II] images are shown in Figure 1. The images contain three 
very bright Milky Way stars directly superposed on the disk. We masked them out in IRAC 
and 24 yum images, as well as additional point sources as described in Section 3.1 below 
(Figure 4b). Foreground stars will be removed more rigorously in a later paper. Fainter, 
indistinguishable stars are unhkely to affect significantly the results presented here. 

3. Physical Content in Mid-Infrcired Images 

Whereas at A > 24 fim the emission is dominated by interstellar dust, the mid-infrared 
marks the transition from stellar to interstellar emission. The appearance of the galaxy shifts 
from resembling the stellar disk in the 3.6 /^m image to an ISM-like appearance at 8 /xm. We 
assume conservatively the 3.6 /im emission is all due to stars, and extrapolate this component 
to longer wavelengths using stellar population modeling from Starburst99 (Leitherer et al. 
1999); the stellar Spectral Energy Distribution (SED) longward of 2 //m depends very httle on 
star formation history or metaUicity, as illustrated in Figure 2. The decomposition amounts 
to scaling the 3.6 fim map by 0.596, 0.399, 0.232 and 0.032 respectively for 4.5, 5.8, 8 and 
24 yum, and subtracting these from the observed maps pixel by pixel to yield what we will 
refer to as the "dust maps" shown in Figure 1. Based on this decomposition, the 24 fim map 
is probably > 98% interstellar emission globally, and > 93% interstellar in the inner arcmin, 
where the stellar contribution is the greatest. The 8 fim map is > 81% interstellar emission. 

This method potentially overestimates the dust component by underestimating the red- 
der SED of very low mass stars. However the model reproduces adequately the 4.5 fim/3.6 fim 
color of the nuclear stellar cluster, of the emission averaged in the inner arcmin (Figure 2), 
and of the diffuse emission outside several disk H II regions. 

On the other hand, the method probably underestimates the dust component by as- 
suming that the 3.6 fim map contains only photospheric emission, ignoring evidence that the 
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hot dust component stretches down to A < 3 //m (Bernard et al. 1994; Lu et al. 2003; Hunt 
et al. 2002). From the ISO data in Figure 6 of Lu et ah (2003), this underestimate may 
reach a factor of up to 2.6, 10-20% and 3-4%, respectively at 4.5, 5.8 and 8 fj,m. In view of 
these opposing potential biases, our conservative approach should be reliable qualitatively 
for NGC 300, in spite of localized departures due to variations in stellar population colors 
or dust extinction. 



3.1. Interstellar Dust Emission 

In order to constrain the origin of the 4.5/im non-stellar emission, we identified pix- 
els where the excess is above the 3a significance threshold in the 4.5/im "dust map", and 
examined their colors and locations. In scatter plots of these pixel values vs. Ha or vs. 
8 yUm dust emission, a correlation branch is evident in both cases with moderate slopes in 
fu{4:.5 iim)/ f^{8 iim) and /;y(4.5 ixm)/}la (Figure 3), as well as a branch with very steep, al- 
most vertical slopes. The pixels in the latter branch are strongly associated with the 3.6 //m 
emission, including the brightest foreground star and a large number of point sources spread 
almost uniformly accross the whole map. We use this excess map as a mask to efficiently 
remove foreground stars, with additional hand masking of residues around very bright stars. 

On the other hand, the branches with moderate 4.5//m/8//m and 4.5 //m/Ha ratios 
map into bright star formation regions, as might be expected (Figure 4a). The global ratio 
vfyiyA.b /im)/z//,^(8 /xm) for the selected dust emission amounts to ~ 4%. However, summing 
the pixels in star forming regions with large 4.5 /im excess, dust alone accounts for 17% 
of the total 4.5 /im emission, probably more if dust contribution to the 3.6 /im emission 
were accounted for. Two uncertainties remain which will be addressed in future papers: 
Whether redder SEDs of young stars affect significantly the results in star forming regions; 
and whether the 4.5/xm dust associates more closely with ionized regions or with Photo- 
Dissociation Regions (PDR). 

The dust maps at 5.8 and 8 yum should be dominated by the Aromatic Features in 
Emission (AFE) proposed by Puget & Leger (1989)) to explain IRAS data, and well studied 
with ISO (Boulanger et al. 1998; Helou et al. 2000; Lu et al. 2003). The value and variability 
of the ratios between these bands are of considerable interest in constraining the physical 
state of the emitters (e.g. de Frees et al. 1993). Lu et al. (2003) found the AFE proffies 
among galaxies constant to within measurement uncertainty of 15 to 25% in the range 5 to 
12/xm; the Spitzer images can extend this finding to much greater accuracy, and could point 
to variations within disks on a variety of scales. 
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We find that ufi,{h.'&^m)/vfi,{^^m) for dust emission alone is constant to within the 
uncertainties, with a value of 0.50 ±0.09 (3(t) over most of the area where surface brightness 
is above the Gcr level in each of the 5.8 /xm and 8 //m dust maps. 

3.2. Comparison to Ha 

Figure 4c shows a comparison between Ha and 8 jjm maps at a resolution of 2 x FWHM ~ 
5.15" ~ 50 pc. In the outer disk, several H II regions as well as more complex structures are 
visible in the Ha map, with the corresponding 8 /im profile flatter and sometimes significantly 
offset from the Ha peak. In spite of the complex structure. Figure 4c-d clearly shows the 8 /xm 
emission highlighting the rims of H II regions. The ratio 8 /im/Ha goes up by a factor five 
to ten from inside to just outside the H II region. This ratio may increase by another factor 
of 2-4 in diffuse regions, but this is uncertain because of limited sensitivity and inaccurate 
continuum subtraction in the Ha map. While this behavior is not a complete surprise, the 
Spitzer data are the first to show it so clearly in a nearby galaxy. Boulanger et al. (1990), 
Sellgren et al. (1990), Giard et al. (1994) and others have reported sharp boundaries to the 
AFE, with peaks on molecular cloud surfaces outside ionized regions, and related this to a 
combination of UV excitation and photo-processing of the AFE carriers. Helou et al. (2001) 
also speculated that AFE arise primarily in PDRs based on physical arguments rather than 
direct observations. Inside the H II regions, the Aromatics are depressed probably because 
of destruction by the ionizing UV (Boulanger et al. 1990). 

The sharp contrasts near H II regions between 8 ^m. and Ha are only a small-scale 
property. Summing over scales larger than about 1 kpc, the ratio between the two tracers is 
constant to about 40%. The global ratio of 8 /im to unextincted Ha is roughly 30% smaller 
than the average measured in disks of spiral galaxies by Roussel et al. (2001a), which had 
itself a dispersion of 50% among disks. For this comparison, the total Ha-|-[N II] flux from 
Hoopes et al. (1996) was corrected for an extinction 74(Ha) — 0.37 and an Ha/(Ha-|-[N II]) 
ratio of 0.91, averaged from values of individual H II regions in Webster & Smith (1983). 
The 8 jivcv IRAC flux was converted to 7 //m flux in the ISO filter LW2 assuming the same 
spectral shape as in the disk of M83 (Roussel et al. 2001b), which was verified to reproduce 
accurately the 5.8 //m/8/xm dust ratio of NGC 300. 

The star formation rate in the disk of NGC 300 is 0.08-0.11 M^y^-^ from the 8//m 
emission (as calibrated for solar- metallicity spiral disks by Roussel et al. 2001a), and ~ 
0.14 M0yr~^ from Ha (as calibrated by Kennicutt 1998). The agreement points to 8//m 
under these conditions as a viable proxy for ionizing flux, and therefore for on-going massive 
star formation rate. Since the association between 8 fim. and Ha is not at the atomic process 
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level however, it will be subject to geometric and other effects, and the use of 8 //m as such 
a proxy is unlikely to be valid under all conditions. 

Figure 4e-f also shows the 24 /im emission is more peaked than 8 /im in star forming 
regions. The spatial resolution is insufficient to distinguish whether PDRs or ionized regions 
are more closely associated with the 24 /im peaks. This confirms earher suggestions based 
on global SED analysis that 24 /im traces heating by the youngest and most massive stars 
in a galaxy (Helou 2000). 

4. Mid-IR to Far-IR SED 

In order to study the SED at A > 8 /im, we identified nineteen emission regions in 
the dust maps, most centered on a local peak (Figure 51). After smoothing all maps to 

the resolution of 160 fira, we measured the emission of each region in a circular aperture 
2.5 X FWHM(160/im) = 95" ~ 1 kpc in diameter. The resulting photometry is plotted 
in Figure 5 as various color ratios against /j,(8 //m)//j,(24yum), and compared to the ratios 
expected from pre-Spitzer model SEDs of galaxies (Dale & Helou 2002). These arc single- 
parameter models of dust emission integrated over whole galaxies, and incorporate a power- 
law distribution of dust over heating intensities. The data and model agree in the trends, 
though there are systematic offsets in the values. These offsets are most simply understood 
as a systematic discrepancy of ~ 30% at 160 //m and ~ 15% at 70 //m between model and 
measurement, assuming reliable photometry at 8 and 24 /xm. In view of the simplicity of the 
models, photometric uncertainties (especially non-linearity effects) and source extent, the 
agreement between model and data is quite satisfactory, and suggests that beyond a certain 
scale portions of galaxies may approach full galaxies in their behavior. 

The 160 fim. and less so the 70 //m maps display a striking large halo which echoes several 
features of the HI envelope. It is unclear at this point whether this reflects unusually cold 
material in the outer disk, or whether the larger 160 yum beam is simply more sensitive to 
low surface brightness features. 

5. Sumrriciry 

The Spitzer Space Telescope data allow us more clearly than ever to dissect the inter- 
stellar dust emission components in NGC 300 and to relate them to star formation. The 
8 fim AFE highlights the rims of the H II regions defined by Ha, confirming it is more closely 
associated with PDRs than with ionized regions. When averaging over regions larger than 
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~ 1 kpc, the ratio of Ha to AFE emission in NGC 300 is consistent with the disk value 
observed in other galaxies. This confirms AFE as a convenient massive star formation rate 
estimator in disks of galaxies. The mid-IR spectral signature of the AFE is invariant as 
measured by the ratio /,^(5.8/im)//i^(8/im), which is essentially constant at ~ 50 pc reso- 
lution at the ~ 20% {3a) measurement accuracy over most of the disk of NGC 300. The 
24 nm emission is more strongly peaked than 8 /xm in star forming regions. Interstellar dust 
emission at 4.5 //m is confirmed at about the expected amplitude; for the first time, its spa- 
tial distribution associates it clearly with star forming regions. The data cannot distinguish 
between PDRs and H II regions as the primary origin for either 4.5 or 24 /im dust emission. 
The SED models of (Dale Sz Helou 2002) apply to ~kpc-sized portions in the disk to within 
the current photometry and calibration uncertainties of 20-30%. 

The Spitzer Space Telescope is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Adminis- 
tration. This research has made use of the NASA/IPAC Extragalactic Database which is 
operated by JPL/Caltech, under contract with NASA. 
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Fig. 1. — Images of NGC 300 obtained with the Spitzer Space Telescope and other sources. 
The left-most column shows observed maps at (a) 3.6 /im, (b) 4.5 /im, (c) 5.8 /xm, (d) 8/im. 
The center column shows (e) BVRHa composite image (credit MPG/ESO), where B is 
coded in blue, V in green, R and Ha in red; (f) Ha map described in Section 2; (g) dust 
map at 5.8 ^m; (h) dust map at 8 /im. The right-most column shows (i) H I map from Puche 
et al. (1990); and observed maps at (j) 24 /im; (k) 70 /xm; (1) 160 /xm, with photometric 
apertures used in Figure 5 superposed. Each frame contains the peak surface brightness 
in the background-subtracted image; all images are in logarithmic scale, except the 160 /xm 
image which is linear. North is about 40° clockwise from vertical up, and each frame is 
approximately 20' on a side. 
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Fig. 2. — SED decomposition: The SEDs of NGC 300 in three different apertures are shown, 
with all fluxes normalized to 3.6 /xm: the nuclear stellar cluster (clearly contaminated by 
dust emission judging from morphology in the IRAC maps), the inner arcmin, where dust 
contribution is relatively small, and the whole galaxy. The thick dashed lines trace model 
photospheric SEDs from 2 to 24 /xm, derived from Starburst99 for a range of star formation 
histories and two metallicities, and transformed to IRAC filter fluxes. The horizontal error 
bars represent the filter widths at half maximum transmission. The Ks fiuxes are lower 
hmits, because the 2MASS map is not sensitive to the diffuse emission. 
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Fig. 3. — Scatter plot of excess /,y(4.5/im) versus dust-only f^{8fim), normalized by the 
noise in the original maps. Individual pixels are plotted, after convolution of all maps to a 

common angular resolution. The correlation branch was isolated by selecting pixels below 
the solid line (corresponding to Figure 4a). The uncorrclated emission, extending to much 
higher brightnesses, was cut here for clarity. The scatter plot versus Ha, not included, shows 
a similar correlation branch, most pixels of which map into the 8 /.tm correlation branch. 




Fig. 4. — Top: Distribution of excess 4.5/im emission, /^(4.5/im), in units of the noise in 
the initial 4.5/im map, (T4.5. (a): Part of /^(4.5/im) which is correlated with 8 /im and Ha 
(950 pixels G [3; 19] (T4.5. (b): Part of /^(4.5 /im) which is not correlated with 8 fim nor Ha, 
in logarithmic scale (4583 pixels G [3; 8150] (T4.5). The contours are Ha isophotes (a) and 
the 4cr isophote of the 8 fim dust map (b) . Bottom: Comparison of Ha and dust emission: 
(c) log(/y(8/im) / Ha) (95% of pixels span 1.2 dex), showing low ratios inside prominent 
H II regions; (e) log(/y(8/im) / /y(24/im)) (95% of pixels span the linear range [0.6; 2.2]), 
showing again the distinctly low color ratios of H II regions; (d) 8 /xm (red) and Ha (green) 
composite image; (f) 8 /im (red) and 24 /im (green) composite image. 
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Fig. 5. — Color-color plot for all Spitzer bands at A > 8 yum, for the regions measured with an 
aperture corresponding to ~ 1 kpc in the plane of NGC 300, with the Dale & Helou (2002) 
model predictions superposed as curves. Each curve has one symbol at x=2.28 to indicate 
the corresponding color ratio. We corrected the 8 fim fluxes by a factor of a/0.69 since the 
aperture is intermediate in size between the nominal point-source aperture and infinitely 
extended scales. The one point with error bars illustrates the uncertainties of 25% on the x- 
axis and 35% on the y-axis. Squares are for 160 /im/8 fim; X signs are for 160 /im/24 /^m; plus 
signs are for 70 /im/8/im; triangles are for 70 /im/24 fim; and circles are for 160 /im/70 /im. 
The systematic offset at the longer wavelength bands is discussed in the text. 



